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The effects of a coupling between tfie quantized mecfianical vibrations of a quantum dot and 
coherent tunnehng of electrons through a single level in the dot are studied. The equation of motion 
for the reduced density operator describing the vibrational degree of freedom is obtained. It is 
shown that that the expectation value of the displacement is an oscillating function of time with an 
exponentially increasing amplitude, which is the signature of a quantum shuttle instability. 
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Nanoelectromechanical systems (NEMS), where elec- 
tronic and mechanical degrees of freedom are couple] 
have been attracting a great deal of attention recentlyol 
An important example of such a system is the nanoelec- 
tromechanical single-electron transistor (NEM-SET) ~ 
a structure where the movable conducting island is elas- 
tically coupled to the electrodes (Fig.l). 
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FIG. 1: Model shuttle system consisting of a movable con- 
ducting island placed between two leads. An effective elastic 
force acting on the dot from the leads is described by the 
parabolic potential. 



In Ref. I it was shown, that the NEM-SET becomes 
unstable with respect to the development of periodic me- 
chanical motion if a large enough bias voltage is applied 
between the leads. This phenomenon is usually referred 
to as a shuttle instability. The key issue in Ref. ^ was 
that the charge on the island, q{t), is correlated with 
its velocity, x{t), in such a way that the time average 
q{t)x{t) > even if i{t) = (see the review Ref. |). A 
classical theory of the shuttle instabilityH was based on 
the assumptions that both the charge on the island and 
its trajectory are well defined quantities. 

A decrease of the island size should result in modifi- 
cations of the electromechanical phenomena in a NEM- 
SET as more quantum mechanical effects come into play. 
There are two different types of quantum effects which 
manifest itself as the island size decreases. The first one 
is the discreteness of the energy spectrum. The electron 
energy level spacing in a nanometer-size grain is of the 
order of 10 K and resonant tunneling effects become es- 
sential at small enough temperatures. In this case the 



characteristic de Broglie wave length associated with the 
island can still be much shorter than the length scale of 
the spatial variations of the "mechanical" potential. If so, 
the motion of the island can be treated classically. The 
NEM-SET in this regime has been studied theoretically 
in Ref. and the conditions for the shuttle instability to 
appear have been found. 

Diminishing the size of the island further results in the 
quantization of the mechanical motion of the island. As 
a result not only the charge on the island but also its tra- 
jectory experience strong quantum fluctuations and the 
picture of the shuttle instability, which was elaborated in 
Refs. H and |7| is no longer valid. 

A NEM-SET system in the regime of quantized me- 
chanical motion of the central island was studied theoret- 
ically in Ref. ^. It was assumed that the phase breaking 
processes are strong enough to make the density matrix 
diagonal in the representation of the eigenstates of the 
quantum oscillator Hamiltonian which describes the 
mechanical subsystem. At the same time it is well known 
that the expectation value of the displacement operator 
in the eigenstates of a quantum oscillator are zero while it 
is the coherent state, which is a coherent superposition of 
all |n), that approaches a description of the classical pe- 
riodic motion of the oscillator as ?i ^ 0. The interesting 
question arises, therefore, how the expectation value of 
displacement operator evolves in time and what state re- 
sults when the formal condition of the shuttle instability 
are satisfied for the quantum NEM-SET. 

In this article we will show that if a large enough bias 
voltage is applied between the leads, then the quantum 
state of the central island of the NEM-SET evolves in 
such a way that the expectation value of the displacement 
operator is not identically zero. Rather, it oscillates in 
time with an increasing amplitudeO. This results shows 
that the shuttle instability is a fundamental phenomenon 
which exists even when the trajectory of the island and 
the charge on it are no longer well-defined. 

We use the following Hamiltonian to model our system 

-ff = ^ eaka'lk^ak + [eo - e£X] c^c + 
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Y,T^{X){ai,c + c^a^k)- (1) 



The first term in the Hamiltonian describes the electrons 
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in the electrodes, the second term relates to the single en- 
ergy level in the central island, the third and forth terms 
to the quantized vibrational degree of freedom associated 
with center-of-mass motion of the central island and the 
last term describes tunneling between the electrodes and 
the island. All energies are measured from the Fermi en- 
ergy of the leads. Here we assume that only one single 
electron state is available in the central island and that 
the electrons in each electrode are non-interacting with 
a constant density of states. 

Let us introduce dimensionless operators for displace- 
ment, X = X/rQ, and momentum, p = r^P/h, where 
ro = \/Ti/ {Mwq), and then measure all lengths in units 
of ro and all energies in units of Tiwq . 

In order to to transfer the x-dependence of the island 
energy level to the term describing tunneling, we make 
the unitary transformation 



H = UHU 



U 



The Hamiltonian can now be written as 



(2) 



(3) 



where 



= ^ ^akO^ak^-ak + C. = Ha + , (4) 

a, A: 



n EE J2To.{nial,c + h.c.] , 



(6) 



with ria = e"^PTa{x), io = €a - MwlcP /2 and d = 
e£/{Mwl). 

The evolution of the system is described by the 
Liouville-von Neumann equation for the total density op- 
erator (t: 



idtcj = [H, a] = [He + Hy, a] + [n, a] 



(7) 



It is convenient to transform Eq. into the interac- 
tion picture, defined by 



Then 



where 



^tiH,+H^)t 



n{t)e 



(8) 



(9) 
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In what follows we will assume that electrons in the 
leads are weakly coupled to the rest of the system and 
that the leads are so large that their statistical properties 
are unaffected by the weak coupling. Then the following 
approximation can be made 



a(t) W p(t) ® CTL (g) cr_R . 



(11) 



Since we are interested in the evolution of the variables 
describing the central island we need an equation of mo- 
tion only for the reduced density operator. If we trace 
over the electrodes degrees of freedom directly in Eq. (H) 
we obtain zero in the RHS, which means that the effect 
is of the higher order with respect to Vl. 

If we formally integrate both sides of Eq. (M) and sub- 
stitute the result back into the RHS of Eq. we obtain 
an integro-differential equation with RHS of the second 
order with respect to fi. 



dta[t) = -i[h{t),5{{))] 



dti [n{t), [n{ti),a{ti 



(12) 

We trace out the electrode degrees of freedom and get 
the integro-differential equation for the reduced density 
operator 



dtp{t) = -Tl-a 



dti 



(13) 



where p{t) = Tia {^{t)} is the reduced density operator 
in the interaction representation and the trace is over the 
electrode degrees of freedom. The first term in Eq. ( p^ ) 
gives zero because we choose such initial condition that 

cr(0) = p(0) (g) (TL ® (TR. 

By using Eq. (Ww we obtain 



dtPt 



E 



de 2?n 



dti e 



ie(t-ti) 



+ [pt,{n^~c'^)t,{nli)t - {nU)Mn^~c^)t,] f~] 

+h.c., (14) 

where /+ = /„ = [1 + e'^^'--^^^]-\ / - = 1 - /„ and 23, 
is the density of states in the corresponding lead. 

This integro-differential equation becomes a differen- 
tial equation if we consider the case of zero temperature, 
T ~ Q, and large bias voltage, eV oo, 

dtpt = 7:VL [2inLC^)tpt{nld)t - [inld)t{nL£^)t,pt} 



To describe the evolution of the oscillator variables we 
need only /5o,o =< 0|p|0 > and pi^i =< l|j5|l >, where 
1 >= ct|0 >. The corresponding equations for /5o,o and 
pi^i are given by 

dtPo,o = -ttPl |fi[fiL,po,o} + 27r2?;j fijjpi4fifl,(16) 

dtpi,i = -7rPi{0;jfiJ^,pi4} +27r23if]ipo.of^I(17) 

It is convenient to change from po.o and pi^i to Rq and 
Ri given by 



i?o(i) 
Ri{t) 



= e ^''e po,oe " e '^'^ , 



(18) 
(19) 
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Then 

dtRo = -i 

dfRi = — i 



a2 



i + 



^'2 



i{fi(a;),i?o} 



,Ri 



Tr{x)Ri^^Tr{x), (20) 

-i{ffl(x),i?i} 



+ ^TL{x)Ro^JrLix), (21) 

where f L,fl,(a;) = TLMx+d/2) and rc«(x) = 2TrVaT^ix). 
These two equations completely describe the evolution of 
the vibrational degree of freedom. 

By using Eqs. (|2^) and (^ij) we can obtain the equa- 
tions of motion for any momenta with respect to the den- 
sity operators R+ = Rq + Ri and i?_ = Rq — Ri. We 
expand Tl^rIx) to first order with respect to the dis- 
placement, Tl,r{x) — {1^x/X)Tl,r{0), where A is char- 
acteristic tunneling length, and leave only terms of the 
first order with respect to A~^. 

In the quasisymmetric case, ^^{0) — Tji{0), the equa- 
tions for n_, x^ and are decoupled from the rest, 



x+ = P+, 

p+ = —x+ — dn- . 

2V 



A 



■x+ 



(22) 
(23) 

(24) 



and equations for the second momenta are decoupled 
from the higher momenta, 



dt {{p,x})_ 



-{{p,x})^~2dp^, 
-2{x^)^ + 2{p^)^-2dx. 



p_ — Tx^ 



2T 
T 



{x-)_ 



(25) 
(26) 
(27) 

(28) 



P- = -a;_-d-fp_ + -({p,a;})+ ,(29) 

where (•)_(_ = Tr{i?±»}, n_ = x± = {x)j_ and 

p± = (p)^ andf = fi(0)+fK(0). 

The characteristic equation for the system of Eqs. ( |2^ , 
P^,p3) is given by 



(a2-t-l)(a-|-f) + 27 = 0, 7 = < 1 



A 



and has three roots 



ai 



27 



r2 + i 



(30) 



(31) 



l-iT 



as 



(32) 



The first root is a negative real number, which corre- 
sponds to a solution which exponentially goes to zero. 
The last two roots have positive real parts and non-zero 
imaginary parts, which gives rise to oscillating solutions 
with exponentially increasing amplitudes. The expecta- 
tion value x{t) = Tr {U^(j{t)Ux} of the displacement op- 
erator X depends on x^ as follows: x{t) = x+ +d/2. This 
means that when we apply a high enough bias voltage be- 
tween the leads the expectation value of the displacement 
starts to oscillate with an increasing amplitude with re- 
spect to the point x = d/2 of the original coordinate 
system. 

The characteristic equation for the system of Eqs. ( |2^ , 
||,|3,|8|,|9l) is given by 

a[A + a^] [1 -I- (f -I- af] + 2j[Afa + 5a^ - 4] ^ , (33) 



and has five roots 



ai 



27 



a2 
a4 



r2 

2i 4 



1 



27 



as 



l-iT 
— r + i , as « a! 



(34) 

(35) 
(36) 



The first root gives a non-oscillatory solution of expo- 
nentially growing amplitude. The remaining four roots 
correspond to oscillatory solutions with increasing {a2 
and aa) and decreasing amplitudes (04 and as), respec- 
tively. Thus the energy of the oscillator, which is given 
by the sum of the two second momenta {x^) ^ and (p^)^, 
exponentially grows with time. 

The importance of the position dependence of Ta{x) 
can be seen if we let A ^ 00 in the above treatment. 
Then we find that the energy grows linearly with time, 
while the average displacement does not grow. 

In conclusion, we have studied a quantum shuttle sys- 
tem where the quantized mechanical vibrations of a quan- 
tum dot is coupled to coherent tunneling of electrons 
through a single level in the dot. We have obtained the 
equation of motion for the reduced density operator de- 
scribing the vibrational degree of freedom. It was shown 
that the expectation value of the displacement is an oscil- 
lating function of time with an exponentially increasing 
amplitude, which is the signature of a quantum shuttle 
instability. 
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